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The behavior of dc subgap current in magnetic quantum point contact is discussed for the case 
of low-transparency junction with different tunnel probabilities for spin-up (-D|) and spin-down 
electrons. Due to the presence of Andreev bound states ±eo in the system the positions of 
subgap electric current steps eVn = (A ± eo)/n are split at temperature T =^ with respect to the 
nonmagnetic result eVn = 2A/n. It is found that under the condition _D| ^ the spin current also 
manifests subgap structure, but only for odd values of n. The split steps corresponding to n = 1,2 
in subgap electric and spin currents are analytically calculated and the following steps are described 
qualitatively. 
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It is well established now that subharmonic gap struc- 
ture in superconducting weak links can be described in 
terms of multiple Andreev reflection (MAR). The concept 
of MAR was first introduced^ for SNS (superconductor 
- normal metal - superconductor) junctions and then ex- 
tended to include the effect of resistance of SN interface^ 
(OTBK theory). Then it was shown^ that this mech- 
anism is the reason for subharmonic gap structures in 
junctions of any type, and specifically in superconduct- 
ing point contacts. As long as the junction is short on the 
scale of the coherence length the microscopic details of 
the junction are irrelevant and for the nonmagnetic junc- 
tion the only parameter is the transmission probability 
of the barrier D. An exact theory of MAR for short non- 
magnetic junctions between clean superconductors has 
been developed in Refs00 on the basis of scattering am- 
plitudes method and in Ref0 making use of nonequilib- 
rium Green function techniques. It was obtained that 
for this case MAR manifests itself in dc current as cur- 
rent steps at voltages eVn — 2A/n, n = 1,2,.... The 
magnitude of the current onset at Vn is proportional to 
Z3" . The theoretical results are consistent with the exper- 
imental situation in atomic-size superconducting point 
contact a''i^i^i^° . The subharmonic gap structure due to 
photon-assisted MAR in the presence of microwave radi- 
ation has also been studied and theoretical predictions^ 
are confirmed by the experimental observationsi^. 

Arnold^ was the first to take into account the phase 
coherence of subsequent Andreev reflections, which were 
absent in the OTBK theory. The quantum coherence in 
short ballistic junctions lead to ac Josephson effec t^i^i^ . 
In addition, in long ballistic SNS junctions coherence ef- 
fects give rise to resonant structures in the dc current 
due to Andreev quantization^^. 

MAR in diffusive nonmagnetic junctions has also been 
investigated in details. Short diffusive junctions with 
Thouless energy Exh ^ A can be described by coher- 
ent MAR theory and show subharmonic gap structures 
at eVn = 2A/nM^. In long diffusive junctions with 
Exh ^ A additional peaks in the differential conductance 
at eVn ~ 2(A ± ETh)/2n are predictedi^. This result is 



in qualitative agreement with the experimental datciii. 
The incoherent MAR regime with eV, A ^ Eth has also 
been studiediSiiS and it was found that subharmonic gap 
structure shows qualitatively different behavior for even 
and odd subharmonic gap structures. 

Magnetic interfaces can not be characterized by 
the only parameter D and Andreev bound states are 
known to take place on the magnetic surface of s-wave 
superconductor—. So it is natural that new qualitative 
features with respect to nonmagnetic case appear in dc 
current even in the simplest models of a short magnetic 
junctionSiiSS. In Refl2a the dc-current-voltage charac- 
teristic of magnetic quantum point contact between two 
superconductors is studied theoretically for the model 
where magnetic quantum dot is characterized by two pa- 
rameters: transparency D and the spin-mixing angle 0. 
It is found that subgap current steps take place at volt- 
ages eVn = A(l -I- cos 6/2)/n. 

In this paper the dc subgap current due to MAR pro- 
cesses in short ballistic quantum point contact for the 
most general case of symmetric magnetic interface includ- 
ing different transparencies for spin-up and spin-down 
quasiparticles (£>! 7^ ^i) is studied theoretically. I con- 
sider the case of small and of the same order of 
magnitude and find that at T 7^ current steps occur 
at voltages eVn = A(l ± cos8/2)/n. This splitting dis- 
appears in the limit T — > and the result obtained in 
Ref221for the current step positions is recovered. Also 
I found that the difference between D-f and is very 
important. It leads to non-zero spin current which has 
subharmonic gap structure different from that one for the 
dc electric current. 

The theoretical analysis is based on the non- 
equilibrium quasiclassical theory of superconductivity 
in terms of Riccati amplitudes developed by EschrigS^ 
and genera lized for the case of magnetic interfaces in 
Refs l20l23 . Now I briefly outline this formalism. The 
fundamental quantity in non-equilibrium quasiclassical 
theory of superconductivity is the quasiclassical Green's 
function g — g{pf,R,e,t). It is a 8 x 8 matrix form 
in the product space of Keldysh, particle-hole and spin 
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variables. In general the quasiclassical Green's functions 
are depend on space i2, time t variables, the direction of 
quasiparticle Fermi momentum pf and the excitation en- 
ergy e. In our case of one- mode quantum point contact 
the problem is effectively one-dimensional and R = x, 
where x - is the coordinate measured along the normal 
to the junction. The momentum pf has only two values, 
which correspond to incoming and outgoing trajectories. 

The electric and spin currents should be calculated via 
Keldysh part of the quasiclassical Green's function. For 
the one-mode quantum point contact the corresponding 
expression for the charge current reads as follows 

+00 
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t'so-o (5^(P/,a;,e,i) - g^(p^,x,e,i)) , (1) 



where e is the electron charge and h ~ I throughout the 



paper. g^{pf,x, e, is a 4 x 4 Keldysh Green's function 
in the product space of particle- hole and spin variables. 
Pf stands for incoming quasiparticle trajectories and p^ 
for the outgoing ones, fj and ai are Pauli matrices in 
particle-hole and spin spaces, respectively. The spin cur- 
rent /s'^ can be calculated making use of Eq.(Q with 
the substitution (T3 for (Tq. s*^ = 1/2 is the electron spin. 

Quasiclassical Green's function g can be expressed in 
terms of Riccati coherence functions ■y^'^ and ^^--^^ 
which measure the relative amplitudes for normal-state 
quasiparticle and quasihole excitations, and distribution 
functions and x^ . All these functions are 2x2 ma- 
trices in spin space and depend on (py , x, e, t). 

In this paper I consider a short junction with the char- 
acteristic size d <C ^0- But it is assumed that despite the 
small size of the grain charging effects can be neglected. 
For definiteness the current is calculated on the left side 
of the interface. Keldysh Green's function for incoming 
trajectory is parameterized by 
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Here subscript 1 means that the corresponding functions 
should be taken at a; = —0, argument p/ of all the Riccati 
functions is omitted for brevity. The product (g) of two 
functions of energy and time is defined by the noncommu- 
tative convolution A(S)B = e'^^^^^ '^^^^^ A{e, t)B{e, t). 
Keldysh Green's function gi{p^) for the outgoing tra- 
jectory can be obtained from Eqs. |(5J, ©by the 
substitution (7f,7f,if)(P/) -> (f f , f f , If )(£p and 
(ff,fii,lf)(p/)-.(^f,7^5f)(£p. 

Riccati coherence and distribution functions obey 
Riccati- type transport equations. The equations are 
given m RefsIHH so I do not write them here. The 
quantities (7^ 7^ if , if ), denoted by lower case 
symbols, are obtained by solving the Riccati equations 
for the appropriate trajectory with the asymptotic con- 
ditions, which are for spin-singlet s-wave superconductor 
as follows 

(A -2ieVi rt 
|e| < A 



xt{e) = {l-\tr({e-eVi,r),t)\'')i^n\,'-^^ , (5) 

where the subscript I, r denotes that the appropriate Ric- 
cati function belongs to the bulk of the left (right) super- 
conductor. A is the bulk absolute value of supercon- 
ducting order parameter for a given temperature, which 
is assumed to be the same in the both superconductors. 
Vi^r is the electric potential in the bulk of left (right) 
superconductor, so V — Vr — Vi is the voltage bias ap- 
plied to the junction. Quantities 7/^^"* and if^ are ob- 
tained from Eqs. I@J, Q respectively by the operation 
d(e, t) ~ a(— e, t)* . 

For analytical consideration superconducting order pa- 
rameter and electric potential are assumed to be spa- 
tially constant in the superconductors. Under this as- 
sumption the voltage drop only occurs at the junc- 
tion region. These simplifications seem to be rea- 
sonable for quantum point contact. Under the as- 
sumptions above the solutions of Riccati equations for 
(^^R^A(^x),^^'^{x),x^ {x),x^ [x)) do not depend on the 
space variable and coincide with the asymptotic condi- 
tions in the corresponding superconductor. 
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The quantities (r^'^^jF^ 
noted by upper case symbols, 

^Ft„A ~R,A 



de- 



ill ' 7i 



are expressed via 
) = (7i^"^.7f"*,^f ) and the 



elements of the interface scattering matrix 5 for the 
normal-state electrons and holes with the energies at 
the Fermi surface. The interface iS-matrix is a unitary 
8x8 matrix in the combined spin, particle-hole and 
directional spaces. The explicit structure of 5-matrix in 
directional space is 



S = 



Sn Si2 

5*21 5; 



22 



(6) 



where matrix Su contains spin-dependent reflection am- 
plitudes of normal-state quasiparticles from the inter- 
face in «-th half-space, while Sij with i ^ j incorporates 
spin-dependent transmission amplitudes of normal-state 
quasiparticles from side j. Each element Sij is a diagonal 

matrix in particle- hole space Sij = Sij{l + T3)/2 + Sij{l — 
T'i)/2. The most general form of S-matrix for a symmet- 
ric magnetic interface without spin-orbit interaction can 
be written asS^: 



Sii — S22 — 



S2 




ie/2 







,-40/2 



ie/2 



-ie/2 



(7) 



(8) 



where + -C^T-i = 1, a = ±1 depending on the par- 
ticular modeP^. Sij ~ Sij in the considered problem. 
The particular expressions for (f f f f Xf- , ) in 

terms of (7^^'"^, if^''^, xf' , if) and 5-matrix elements can 
be found in Ref 24, so I do not write them explicitly. 



Now we can proceed with the current. Substituting 
Keldysh Green's function Q into Eq.Q} it can be seen 
that the current is expressed as a sum over harmonics 
j(^) = 'Y^ Jn&'^^^'^^^ ■ I focus on the dc component and 

N 

only calculate the term corresponding to TV = 0. I-V 
characteristics in charge and spin currents exhibit steps. 
In the tunnel limit one can extract the elementary pro- 
cesses giving rise to each of these steps. In this paper the 
positions, height and shape of several first steps in sub- 
gap I-V characteristic are calculated analytically. The 
linear in transparency thermally-activated term in the 
dc current, which is not related to the steps, has also 
been calculated, but it is not pronounced in the subgap 
region for the considered problem and is not discussed 
below. The part of the dc current corresponding to the 
steps takes the form: 
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A cose/2. 7(e) = 
> A the square root 
by — zsgnevc— A2. 



■ch{sp) 
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(e — i\/A2 — e2)/A, where for |e| 
a/A2 — e2 should be substituted 

The behavior of j'^'* and j'^^ described by Eq. ^ is shown 
in FigUlfor e = 7r/2. 

It is seen from FigHlthat at zero temperature there 
is a shift of well-known step positions in electric current 
for a non-magnetic interlayer: eVn = 2A/n — s- eV^ = 
(A + en)/n. This result is coincide with that one found in 
Refl22. However, at non-zero temperature non-magnetic 
step positions are split: eVn — 2A/n eV^f = (A ± 
£[))/n. It is easily obtained that the additional steps eV~ 
disappear rapidly (cx e~^°/^) when the temperature goes 
down. In the first order of D the analytical expression for 
the electric current similar to Eq. H1U|) has been obtained 



for the case of nonmagnetic junction between two d-wave 
superconductors2&. 

Let us describe all the steps for n > 2 qualitatively. 
The estimate for the magnitude of jn''^*^'' is easily ob- 
tained if one considers physical processes of MAR lead- 
ing to subgap current steps. For the spin-up band the 
bound state energy is £5 = -|-£o and at T = the fol- 
lowing process gives spin- up current at eVj^: a spin-up 
electron tunnels into the bound state from the contin- 
uum states below —A undergoing n — 1 Andreev reflec- 
tions inside the gap. An electron converts into a hole and 
vice versa under each Andreev reflection event. In every 
passing of the electron (hole) through the boundary the 
amplitude of its wave function should be multiplied by 



For the spin-down band with Eb 
the same process takes place with the modification that 
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FIG. 1: (a) The subgap structure of dc electric current (in units of Aoe/fi) is plotted in accordance with Eq.lO for T = Q, O.STc- 
(b) The same for the spin current. 9 = 7r/2, D| = 0.2, = 0.1. Ao = A(T = 0). 



a spin-down electron starts in the bound state and goes 
up into the continuum states at -l-A. So the current 
j\A2k is oc (D^D^f and jTa(2fc+i) « {D^D^fD-^^^. Fi- 
nally, we obtain for j'f!^^^^^ = j„-f ± jni the following re- 
suit: J(tiS (X pT^i)'PT±^i) andj-% « 2pt^i)'> 

■l{2k) — 

For T 7^ there is a non-zero probability to find a spin- 
up electron at Eb = So and a vacant place for a spin-down 
electron at — Eq. So the processes leading to the addi- 
tional steps at eV~ arise: (i) a spin-up electron tunnel- 
ing from = eo to the continuum spectrum at +A and 
(ii) a spin-down electron tunneling from the continuum 
spectrum below —A into the bound state —eg. Conse- 
quently, for T 7^ the splitting of current step positions 
occurs and there are all the steps at eV^ = (A ± eo)/n 
in the charge current and only odd split steps in the spin 
current. 

Of course, the processes of inelastic scattering lead 
to impossibility of observing the subgap structures cor- 
responding to large n due to the loss of coherence if 
the time r — ndef f /vf, which a quasiparticle spends 



in the junction region, exceeds the average time of in- 
elastic scattering Ti„. So the steps with n more than 
Ti„Ao(^o/'^e//) can not be observed. For the case under 
consideration this estimate seems not to be restrictive be- 
cause ^o/deff ~ 1 for a short junction and Ti„Ao due to 
electron-phonon scattering oc {u}d/T)'^{Aq/T) and con- 
siderably exceeds unity for temperatures up to Tc- 

In conclusion, I have presented a theoretical study of 
the dc electric and spin current subgap structure in mag- 
netic point contact. At T 7^ electric current step posi- 
tions eVn = {sq ± A) /n are split compared to the case of 
nonmagnetic junction and at T = they are just shifted: 
eVn = (eo + A)/n. The subgap spin current, taking place 
at D^i Di manifests steps at the same voltages but only 
for odd values of n. The current for steps with n — 1,2 
is calculated analytically and a qualitative description of 
the following steps is presented. 
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